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Abstract 



The competition between antiferromagnetism and the d + id superconducting state 

is studied in a model with near and next near neighbour interactions in the absence of 

a . any on-site repulsion. A mean field study shows that it is possible to have simultaneous 

occurrence of an antiferromagnetic and a singlet d -\- id superconducting state in this 

F^ \ model. In addition, such a coexistence generates a triplet d + id superconducting order 

P3 ' parameter with centre of mass momentum Q = (vr, vr) dynamically having the same 

^ , orbital symmetry as the singlet superconductor. Inclusion of next nearest neighbour 

hopping in the band stabilises the dxy superconducting state away from half filling, the 

topology of the phase diagram, though, remains similar to the near neighbour model. In 

^ . view of the very recent observation of a broad region of coexistence of antiferrmagnetic 

^^O \ and unconventional superconducting states in organic superconductors, the possibility 

zL ' of observation of the triplet state has been outlined. 

O ■ PACS Nos. 74.20.-z, 74.72.-h 
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I. Introduction 



Interests on the interplay between antiferromagnetism and superconductivity date back 

'^ \ quite a while as in certain organic and heavy fermion superconductors, superconductivity is 

known to coexist with an antiferromagnetic (AF) phase at low temperatures |jl|, |^. In a recent 

O ■ study of the organic superconductor fi;-(ET)2Cu[N(CN)2]Cl the most complete phase diagram 

has been obtained as a function of pressure and a region of coexistence of unconventional 

KA \ SC and AF LRO is observed]^]. In the high Tc cuprates, the superconducting state abuts 

J^ \ (albeit with a small gap) the AF state and recent inelastic neutron scattering^, H reveals 

considerable AF fluctuations deep inside the superconducting (SC) state in YBCO. There are, 

actually, quite a few similarities between the organic and high Tc superconductors 0. Indeed, 

antiferromagnetism, in some theories, is considered to lie at the heart of the mechanism that 

drives unconventional superconductivity 0]. 
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Several investigations were carried out^ to study the phase diagram and nature of phase 
transitions between the AF state and the SC state in the context of the organic and heavy 
fermion superconductors. Meintrup et al.[^] proposed a model in the general context of 
antiferromagnetism and superconductivity in which nearest neighbour singlet pairing inter- 
action was shown to accommodate both SC and AF states and their coexistence in certain 
range of parameters. These authors, and several others before them[|iy, |lT], showed that it 
is not necessary to have strong on-site repulsion to generate AF LRO. Correlated models 
with extended range of interactions can produce AF LRO and SC as well. 



In a recently proposed S0(5) theory of superconductivity Zhang |]12| considered a five 
dimensional order parameter space {ipi, ■ ■ ■ ,ip5) with ipi + iip^ being the SC order parameter 
and the remaining three constituting the AF moment. A rotation in this five dimensional or- 
der parameter space, effected by the spin 1, charge 2, vr— operators, (vr^ = J2k 9ik)'^k+Qi'^-ki^ 
where c^(c)'s are the electron creation (destruction) operators and Q = (vr, vr) in d=2, leads 
to the transition from the AF to the SC state. In this theory, the triplet magnetic excitation 
of the quantum disordered phase is identified with this vr-triplet mode in the SC phase. In 



an exact diagonalization study of the t — J model |T3], the dynamical correlation functions of 



the TT— operators have been calculated and found to be non-zero. The existence of both AF 
and triplet pairing amplitude with net momentum Q was reported earlier in the mean-field 
study of a pairing Hamiltonian in the context of the heavy fermions and organic supercon- 



ductors P, |lj], although the conditions under which the triplet amplitude appears and the 
modifications of the phase boundaries due to this triplet amplitude were not dealt with. In 
a recent investigation, Kyung|l5| considered explicit mean-field pairing interactions in the 



singlet and triplet channel with a repulsive on-site interaction to stabilize the d-wave SC 
state (over s-wave) and discussed the coexistence of a dynamically generated triplet SC pair 
amplitude and AF long range order (LRO). In this case, superconductivity is governed by 
the attractive interactions in the appropriate channels while the AF state owes its origin 
primarily to the on-site repulsion in the usual manner. 

We start from a Hamiltonian with nearest and next nearest neighbour interaction and 
consider pairing in the d-wave and d + id-wave in the singlet (and later on in the triplet) 
channel along with the antiferromagnetic LRO. The on-site interaction is assumed to be 
small|]16| (set to zero here) and both the LRO and off diagonal long range order (ODLRO) 
are governed by the same interactions in a manner similar to the case studied by Meintrup 
et. al.0 (where the on-site interaction was absent as well). Using a mean-field analysis, 
these authors studied the coexistence of different singlet SC and AF LRO states in their 



model. Related extended range models have been studied by numerical methods ||T0|] and 
mean-field [|lT| theory earlier, although a detailed study with the possibility of several SC 
symmetries and AF order have not been undertaken. 



There have been suggestions [0 for the existence oi d-\- id state in the high Tc systems 



followed by possible observation in a series of experiments [18, 19|. Recently, it has been 



shown [^] from a renormalization group analysis of the fluctuations that the transition to 



d + id state possesses a stable fixed point. Kino and FukuyamaPT[] considered a model with 
only on-site repulsion for the organic superconductors in the intermediate coupling range 



(typically the on-site interaction is about half the band width). Such a model, although 
accounts for the AF phase and the metal-insulator transition, fails to explain the large 
SC phase observed in^]. Extended range attractive interactions with the right symmetry 
are necessary to obtain these unconventional SC states. Additional processes discussed in 
ref.[|T6|, particularly in the large metallic region above the SC phase, possibly reduce the 
on-site interaction in the organic systems further. 

A preliminary report of the coexistence of and competition between the singlet supercon- 
ducting d^2_y2, dxy, dr^2_y2 + idxy (the so called d + id state) states and AF LRO in a model 
similar to that of Meintrup et. al. with extended range of interaction has been presented 
recently by two of usp2|]. We extend this calculation and show in the present work that in 



the presence of such coexisting singlet SC order parameter and AF LRO, a triplet pairing 
amplitude with centre of mass momentum Q is dynamically generated even if there is no 
explicit interaction in the Hamiltonian in that channel. It is not a-priori obvious that the 
dynamical generation of the triplet amplitude should occur in a model where the AF and 
SC states are governed by combinations of the same interactions. In the model considered 
by Kyung[|I^ the relative strength of these two competing states are governed primarily by 
separate and independent interaction parameters. We also figure out how the phase diagram 
gets modified in the presence of the triplet pairing amplitude in the present model. In section 
II we discuss the model under consideration. Section III concerns with the results, discussion 
and concluding remarks. 

II. Model and Calculations 

The model studied here incorporates antiferrmagnetic LRO and superconductivity and is 
given by the Hamiltonian [^ 

H = ^^kcLck<7+ Yl Vijnianj^„+ ^ VijUiaUja' (1) 

k,<T <ij>c7 «ij»rja' 

where the sum over < ij > extends over near neighbour and << ij » over next near 
neighbour sites. We take Vi and V2 as the corresponding interaction strengths (both Vi and 
V2 are negative) and write ^k = Ck — /W. In the absence of an on-site repulsion, this model is 
perhaps the simplest that produces AFM as well as superconductivity in the d + id channel. 
Elementary physical reasoning shows how an AF state appears in this Hamiltonian. In 
the classical limit (tjj = 0) the Hamiltonian has near-neighbour attractive density-density 
interaction (amongst opposite spins) that leads to AF spin correlation (of Ising symmetry) 
among nearest neighbour spins. An on-site repulsive term would have stabilised this further 
and the region of AF LRO would extend in the phase space. The second neighbour attractive 
density-density correlation term is spin independent and stabilises a d^y order in the quantum 
limit. The regions of stability of dx2_y2, d^y and the d + id state for the range of values of Vi 
and V2 have been discussed in[^]. Extensive literature exists for models in the opposite limit 
of repulsive extended range interactions where the classical limit gives rise to charge density 
waves i23[|. In the absence of V2, Monte-Carlo calculations |T0(] of model (1) shows AF phase 



at all densities. At half-filling the ground state is Neel ordered while away from half-filling 



there is evidence for phase separation between AF ordered and empty domains. Mean-field 
analysis captures much of these features qualitatively, although a realistic description of 
the phase separation eludes such calculations as expected. 

In order to use mean-field description for the symmetry-broken states we define the 



operators corresponding to the singlet and triplet SC order parameters in real space |p!6| 



^i.^ ^ I ^ crCi+5,<^Q _^0(5) and A^^t = jYl Ci+s,aCi-a(f)i5) (2) 

5,(7 5,(T 

where 6 is the usual near-neighbour translation vector and a choice of the form factor 0(5) 
with 0i(5) = 1 for 5 = (±1, 0) and 0i((5) = -1 for 6 = (0, ±1) (with ^S) = for all other 
choice of 6) ensures that the SC OP has d^2_y2 symmetry. For d^y symmetry, one takes the 
form factor (the only non-zero terms) (p as 02(^) = 1 for ^ = (il? =tl) and 02(^) = —1 for 
6 = (±1,=f1). The operator corresponding to the AF order parameter is the well known 
form Z^o-C" 4,0-^4,0- • Writing the AF, singlet and triplet SC order parameters as {a = ±1) 

Y, < (^c{aCi,a >= boe'^-"^ (3a) 

- ^ < aCi+s,aCi-a > (01 ((^) + i(p2iS)) = A, (36) 

^ 5,a 

and 

'^ 5,a 

The SC order parameters are chosen to be of dx2-y2 and d^y symmetries as the interactions 
Vi and V2 are known to favour |]20|, ^ , ^ superconductivity in such orbital symmetries. The 
presence of AF order in the model (1) has already been indicated. The superconducting order 
parameter A^, which is a spin singlet with d-wave orbital symmetry, can have a non-zero 
value when the underlying pairing state is spatially homogeneous. On the other hand, the 
order parameter Af that we write down here should be thought of as a non-zero expectation 
value of the vr— operator of the S0(5) theory and is generated only dynamically EM. It can 



be non-zero provided that the expectation value of the annihilation operators in (3c) changes 
sign on alternating bonds. The resulting Aj is thus a staggered order parameter: it changes 
sign from site to site in the same way as the antiferromagnetic order parameter (3a) does. 

Using the Hartree-Fock approximation with these order parameters and going over to 
Fourier space we get a Hamiltonian in quadratic form as 

Hmf = ^(ek - /i)cLck<7 + &mIZ[(CkTCk+QT - 4iCk+Qi) + h.c] + ^(Z\*(k)c_kiCk| + h.c.) 

k,a k k 

+ IZ[^t(k)(c-kiCk+Qt + c_k-QiCkt) + h.c] (4) 



Here As(k) = ^Ai{coskx — cosky) + i/S.2sinkxSinky and Ai(k) = \/S.ti{coskx — cosky) + 
iAt2sinkxSinky are the pairing amplitudes in the singlet and triplet channel respectively; 
bm = Abo with A = zV\ + ZV2 {z is the coordination number). The tight binding energy 
dispersion on a square lattice that we use involves upto next-near neighbour hopping in 
conformity with the range of interactions considered. The dispersion is ek = —2t{coskx + 
cosky) + At' cosk^cosky where t is the near neighbour and t' is the next near neighbour hopping 
integral. Here fl = ft — nziVi + V2) and we use /x to denote jl in the foregoing analysis. 

The potential and the SC order parameters are expanded p6|, p^ in the usual basis func- 
tions {Bi representation of C4t,) r7i(k) = ^{cosk^— cosky) andr72(k) = sink.j.sinky. Expanding 
in these bases \^(k, k') = X^i ^j^j(k)?7j(k') and A^ = I]k' ^(k, k')rk' = X]i^Ajr7j(k) where 
Aj = Z]k^j(k)rk (and Fk =< Z^o- Cko-c_ko- >)• ^(k, k') comes from the Fourier transform of 
the second and third terms of the Hamiltonian (1) (described in detail in ref. p6[]). 

The order parameters can be treated as variational parameters with the trial Hamiltonian 
Hmf- The corresponding free energy functional is given by 

F = Fo+<H- Hmf >o (5) 

where < ... >o denotes average with respect to po = Gxp[—HMF/kT]/ZQ and Fq = —kTln^Zo). 
The self-consistency equations for the order parameters are obtained by minimising the free- 
energy functional F. 

2iV^^^_ ^ a(k)^E^(k) 2T ^ ' 

1--EE ^-^tanh^^ (7) 

1--EE ^-^tanh^^ (8) 

2N ^ ^ EJ\<:\ IT ^ ' 



2^ ir 7^^- ^7(k) 2T 



A« - -^\Y._iyV.^^,l(^.),AK^^^anK^ 



(9) 



The self-consistency equations (6)-(10) provide a set of five coupled equations to be 
solved numerically. What is interesting is that although we have not included any pairing 
interaction in the triplet channel, the triplet amplitudes Ati and Ai2 are non-zero. The 
simultaneous co-existence of Ai^2 and 60 ensures the existence of this non zero amplitude. 

One could, of course, include additional pairing interactions W\ and W2 explicitly in 
the Hamiltonian (1) in the triplet channels with preferred symmetries and obtain the self- 
consistency equations. (An alternative and quite commonly used approach is to include such 
terms in the mean-field Hamiltonian (4) as in Kato and Machida^] and Kyung|T5|, rather 



than derive them from microscopic interactions). This would modify all the equations (6)- 
(10) but the nature of the phase diagrams remains qualitatively similar (discussed later). For 
the sake of completeness, we write down the equations for the triplet amplitudes including 
Wi and W2 and note that they reduce to (9) and (10) without Wi^2- 

^^^ k 7=+ - I ^^^) 

+V2W2A2A,2vl{k) - f,Abo) + -^^eiW^Anvl{k)'^ ^E^ky^"^^^ ^^^^ 
and 

^*2 = ^E E \w2At2vl{^) + ^^V2A2vl{^){V^W,A,AWM 
2iV ^ ^^+ _ (^ a(k) 

+K2l^2A2A,2r/2^(k) - M&o) + ^^eiW2A,2vi{k)^ E^kJ^"^^^ ^^2) 
The energy eigenvalues E^Ck) used above are 

i^7(k) = [bl + |A.(k)p + |A,(k)p + el + f,' + 27a(k)]V2 (13) 

where a(k) = {{ViWiAiAtivlik) + V2W2A2At2vli^)-f^bmf+elifi^ + \Atik)\^)y/^ with 1^1,2 
set to zero in equations (6)-(10). In order to obtain the phase diagram in the temperature- 
density plane, the particle density n is calculated from n = — ^ as 

n=l + 7^Y: E l^^ + ^,[i-Abo){V^W^A^AWM 

+\/2W-2A2A,2r72(k) - M&o) + 4/"] } ;^^«^^^ (14) 

Note that ^k + ■Ck+Q = when t' = 0. The k-sums run over half the Brillouin zone to 
accommodate the zone folding due to AF state. 

III. Results and Discussion 

It is straightforward to check that these self-consistency equations reduce to simpler and 
well known forms in the limit of pure phases (either AF or SC). Setting Af and A^ zero in 



equation (7), we recover Ey(k) = —ji — 'yJ{b'^ + e^) (where 7, as before, is ±1) and Q;(k) = 
(bm + ^k)- This leads to the well known self consistency equation for AF order parameter 

1/A = — I I]k,7=±i '^ °" Q(k) • Similar reduction occurs in the equations for singlet or triplet 
SC order parameters in the absence of other two. The complete solutions of the non-linear 
coupled set of equations (6)-(ll) have been obtained numerically. Before discussing these 
solutions and the resulting phase diagrams, we examine some of the equations critically. 

6 



The structure of these self-consistency equations for the order parameters lend themselves 
to some interesting conclusions as noted earlier. The amplitude Aji has a finite value even 
when the pairing interaction in the triplet channel with corresponding symmetry is zero, 
provided the AF (bo) and the dx:2_y2 SC order parameters (Ai) are non-zero simultaneously. 
In exactly similar manner At2 gets dynamically generated when both AF and dxy SC order 
parameters (i.e., 60 and A2) appear simultaneously while the pairing interaction(PF2) is zero 
in eqn.(lO). Simultaneous presence of AF and singlet d + id SC state dynamically generates 
the triplet d + id SC state. This is a reflection of the fact that the presence of spin density 
wave order parameter < cj^ |Ck+Q,t > and singlet SC order parameter < cj^ ^cLk,; > can 
lead to a coupling in the triplet channel < cJ^_,_q|C_,^ ^ >. Note that the symmetries of the 
order parameters for singlet and triplet SC states have to be the same. A glance at the 
terms causing the dynamical generation of triplet SC order parameters in eqns. (9) and (10) 
reveals that they contain a factor 77^ (k) of which one 77j(k) term comes from spin singlet 
amplitude and the other comes from the spin triplet term. For the dynamical generation of 
triplet SC order parameters, it is necessary that both of them have the same symmetry or 
at least non-orthogonal. At half filling (/i = 0) these terms responsible for the dynamical 
generation of triplet amplitudes vanish and there will be no triplet SC state in the absence 
of Wi or W2. 

The self-consistency equations for the order parameters are solved numerically for differ- 
ent values of the interaction strengths Vi, V2, Wi and W2 in the presence of nearest {t' = 0) 
and next nearest neighbour hopping (t' 7^ 0). The corresponding phase diagrams are shown 
in Figs. 1-5 in the doping {x = n — 1), temperature (T) plane. In this calculation, all ener- 
gies and temperatures are scaled in units of t. In Fig.l, where Vi = —0.17, V2 = —0.08, 
Wi = W2 = the ground state is antiferromagnetic at half filling for t' = 0. As we move 
slightly away from half filling, a phase appears where the order parameters corresponding to 
AF, dx2_y2-S>C and the vr-triplet SC with dx2_y2 symmetry are simultaneously present. Note 
that the triplet SC phase appears even though the pairing potential in the triplet channel 
(Wi) is zero. This phase is generated dynamically in the presence of the other two phases. 
In the region where different ordered phases coexist, the system actually phase separates [^J: 
there is a first order transition between the SC and the AF states. In the mean-field theory 
there is a single phase boundary that separates the two phases, whereas in an actual sys- 
tem, with long range interactions, there could be multiple domains of one phase in another. 
Finally, far away from half filling we get an SC-only phase having d^i^yi symmetry. 

Fig. 2 describes the phase diagram with a larger value of V2 = —0.32 keeping the other 
parameters same as in Fig.l. The phase diagram has the same topology as in Fig.l, but 
an increased V2 favours the d^y state over the dx2-y2 and hence the three phase region now 
comprises of AF, d^y-'&C and the vr-triplet SC having d^y symmetry . As before, the triplet 
SC phase exists (with a different symmetry compared to Fig.l, being forced by the d^y 
symmetry of the corresponding singlet phase now) even without the pairing interaction W2- 

The generic phase diagram of the model remains similar to either Fig.l or Fig. 2 as the 
ratio y- is changed. As noted earlier |2^, simultaneous appearance of an AFM phase with 



d + id SC can be observed in a narrow regime of parameter space with strongly suppressed 



AFM region. In a similar vein, we observe in Figs. 3-5, a non-zero value of both AFM and 
d + id order parameter (along with the dynamically generated triplet amplitude), when the 
interaction strength corresponding to AFM amplitude is suppressed {A = zVi). The x-T 
phase diagram for Vi = —0.21, V2 = —0.32 with Wi = W2 = and t' = is shown in Fig.3a. 
At half filling the ground state is antiferromagnetic as usual. On increasing the filling slightly 
a phase appears where AF, d^^-y^ + id^y SC and the triplet dx2-y2 + id^y SC amplitudes are 
simultaneously non-zero. The triplet superconducting amplitude is generated dynamically 
and has the same symmetry as of the singlet SC state as expected. On increasing the doping 
we get a singlet dx2_y2 + id^y SC phase. Note that at a higher temperature in the phase 
diagram, there is a region where only the pure dxy phase survives. The same phase diagram 
can be drawn in the temperature-chemical potential plane (Fig. 3b), a situation that obtains 
in some experiments where it is difficult to dope the system while as a function of pressure 
there are interesting phase transitions observed. It is quite interesting to note the similarity 
between the phase diagram shown in Fig. 3b here with Fig.l in Lefebvre et. al.|^. The 
symmetry of the SC phase abutting the AF phase in our model depends on the values of 
the parameters Vi and V2. A dynamical generation of triplet pairing amplitude, therefore, 
remains a distinct possibility in the region of coexistence of AF and SC phases in the organic 
superconductors and further experiments are needed to ascertain this. We did not extend 
our study of the model to the high temperature normal state properties and it remains to 
be observed if strong pairing fluctuations render the single particle properties of that state 
unusual. 

In the phase diagram shown in Fig. 4 pairing potentials for vr triplet SC state are taken 
to be finite. The topology of the phase diagram remains the same as in Fig. 3. The phase 
boundaries shift to provide a larger triplet region only and no separate triplet SC region 
arises even if the values of Wi and W2 are made comparable to those of Vi and V2. 

In the phase diagram of Fig. 5 we introduced a small t' = 0.03 without changing Vi and V2 
and keeping Wi = W2 = 0. Changing the band structure with a non zero t' is known ||25|] to 
favour the dxy state over the d^^^y^ state. But simultaneously, the chemical potential shifts 
for a particular doping on introduction of t'. These two effects act counter to each other 
in the term ^k + ^k+Q (which is — 2/i if t' = 0). As a result, we get a small sliver of pure 
dx2-y2 component before the d + id phase and the coexistence region therefore has AF, Ai 
and the dynamically generated An. At high temperature we get a pure d^y SC phase, below 
which a dx2-y2 + id^y state appears. On increasing the next near neighbour hopping, the dxy 
phase stabilises and the sliver of dx2-y2 state disappears (Fig. 6). There is a combined region 
of AF and d + id state in both singlet and triplet channel. The triplet part is, of course, 
dynamically generated here as Wi = W2 = 0. For the purpose of demonstration, we have 
shrunk the region of pure AF phase around half filling in Fig. 6. We also observe that in our 
model the on-site Coulomb interaction term has been set to zero. Presence of this term, at 
least at the mean-field level, will not change the phase diagram qualitatively[jT5|; although 
the region of stability of the AF phase increases with such a term. 

In conclusion, we have demonstrated the possibility of simultaneous presence of a spin 
density wave and superconductivity in the d + id channel in a model with only extended 



pairing interactions. As it turns out that such simultaneous appearance of AF and singlet SC 
order parameter leads to a spontaneous generation of a triplet SC amplitude with the same 
symmetry even when the corresponding pairing interaction in the triplet channel is absent. 
Though we do not intend to propose the present model for the organic superconductors, 
the nature of phase diagrams we obtained from a generic correlated electronic model that 
produces coexistence of AF and singlet SC state bears similarity to the ones obtained for 
them. Since the dynamical generation of the triplet amplitude rests only on the coexistence of 
AF and singlet SC order, we believe it would be interesting to see if further measurements in 
the organic superconductors like fi;-(ET)2Cu[N(CN)2]Cl reveal the presence of unconventional 
superconductivity in the triplet channel as well. 
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Figure captions 

Fig. 1. Phase diagram in doping x = n — 1 and temperature (T) plane for Vi = —0.17, 
V2 = —0.08, Wi = W2 = and t' = 0. All energies are measured in units of t. 

Fig. 2. Phase diagram in doping and temperature plane for Vi = —0.17, V2 = —0.32, Wi = 
W2 = 0, t' = 0. 

Fig. 3. a) Appearance of the dynamically generated triplet state for Vi = —0.21, V2 = —0.32, 
Wi = W2 = 0, t' = 0. All five amplitudes are non-zero in the region of coexistence. In 
b) is shown the same phase diagram in the T-fi plane. The dashed line represents a 
first order transition, while the solid lines stand for second order transitions. 

Fig. 4. No major change appears in the phase diagram with Wi = —0.21 and W2 = —0.32 in 
comparison to Fig. 3. The triplet phase occupies slightly larger region in the phase 
diagram. 

Fig. 5. A non-zero small t' = 0.03 changes the phase diagram with the sliver of dx2_y2 SC 
state appearing in between with no d^y state in the coexistence region. Both Wi and 
W2 have been kept zero. 

Fig. 6. A larger t' = 0.1 gives rise to the region of coexistence of all five amplitudes again. 
The dxy component stabilises on increasing t'. 
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